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Abstract 

Recently, I. Brevik, J. B. Aarseth, J. S. H0ye, and K. A. Milton [Phys. Rev. E 71, 056101 (2005)] 
adduced arguments against the traditional approach to the thermal Casimir force between real 
metals and in favor of one of the alternative approaches. The latter assume zero contribution from 
the transverse elèctric mode at zero frequency in qualitative disagreement with unity as given by the 
thermal quantum field theory for ideal metals. Those authors claim that their approach is consistent 
with experiments as well as with thermodynamics. We demonstrate that these conclusions are 
incorrect. We show specifically that their results are contradicted by four recent experiments and 
also violate the third law of thermodynamics (the Nernst heat theorem). 

PACS numbers: 42.50.Pq, 03.70.+k, ll.10.Wx, 78.20.Ci 
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The paper is devoted to an important problem which has been actively discussed in the 
literature starting in 2000, and which has created much controversy among various research 
groups. The authors of Ref. try to justify an alternative approach to the theoretical 
description of the thermal Casimir force, which was first suggested in Ref. and later 
supported and further developed in their own Refs. 0,0,0]. (Henceforth, this approach will 
be called BSBAHM after the principal authors Bostròm, Sernelius, Brevik, Aarseth, H0ye, 
and Milton.) However, in Ref. no attention is paid to the serious shortcomings of this 
approach, and the traditional approaches in their discussion are misrepresented. 

The authors of Ref. ÏJ claim that their calculations of the thermal Casimir force between a 
Cu plate and an Au sphere are: a) consistent with current experiments (their wording is "not 
inconsistent"); b) consistent with the third law of thermodynamics. Below we demonstrate 
that both of these conclusions are incorrect. 

a) The first main claim contained in Ref. [l is that the approach by BSBAHM is consistent 
with current experiments (Sec. III). In Ref. 6j the opposite conclusion was drawn, namely 
that this approach is experiment ally excluded. Note that Ref. jfji] contains the results of two 
experiments performed using a micromechanical torsional oscillator: a static measurement 
of the Casimir force between a sphere and a plate, and a dynamic measurement of the 
effective Casimir pressure between two plane parallel plates. The static measurement was 
less precise than the dynamic measurement over a wide separation range. The conclusion 
that the BSBAHM approach is inconsistent with the experimental data in the separation 
region from z = 260 nm to z = 700 nm was made in Ref. 0] on the basis of the dynamic 
measurement. This conclusion is illustrated in Fig. 12 of Ref. where the differences 
between the theoretical P^ 1 ' 1 (calculated in the framework of BSBAHM approach) and 
experimental P^ xp vàlues of the Casimir pressure are plotted as a function of separation for 
one set of measurements containing 235 experimental points. For comparison, in Fig. 11 of 
Ref. | the same differences are plotted for the theoretical Casimir pressures calculated in the 



framework of the impedance approach of Ref. pj . (Recali that this approach, as well as the 
approach using the dielectric function of the plasma model, are called "traditional" since they 
yield results in qualitative agreement with the case of ideal metals; the thermal corrections 
predicted by the alternative approach of Refs. □ DDDD are many times greater at short 
separations.) From these two figures it becomes apparent that the BSBAHM approach is 
excluded by experiment, whereas the impedance approach is consistent with experimental 
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data. 

Surprisingly, in order to demonstrate the consistency of their approach with current ex- 
periments, the authors of Ref. |j] discuss at length only the experiments of Refs. j^, |ç| and 

are at present not the most accurate, and 
were not used in the literature to exclude their approach. As to the dynamic measurement of 
Ref. (which contradicts the BSBAHM approach), only a brief mention of this experiment 
is made at the end of Sec. III, even though this is the most precise current measurement. 
Ref. jl| claims without proof that there are systematic theoretical and experimental uncer- 
tainties in Ref. 6 1 connected, in particular, with a systematic shift of position as discussed 
earlier in Refs. 5, 11 1. It is easy to verify, however, that the reasoning of Refs. ^, 11 1 is 
incapable of avoiding the conclusion of Ref. p that the BSBAHM approach is excluded 



experimentally. According to Refs. |5 , 



ll| , even a small experimental error Az in separation 



distances between two parallel plates (Az = 1 nm in Ref. [6j]) leads to an error in the theo- 
retical Casimir pressure given by AP 1 ^' 1 m —P^^^AAz/z). At short separations this error 



may be rather large. At separations under consideration in Ref. |6] it is, however, much less 
than the discrepàncies between the BSBAHM theory and experiment. Thus, at the shortest 
separation z = 260 nm, AP^ 1 ' 1 = 3.7 mPa, which compares with the 5.5 mPa mean deviation 
of the BSBAHM prediction from experiment as shown in Fig. 12 of Ref. |6j. At separations 
of z = 300, 400, 500 , and 600 nm the above error is equal to 1.5, 0.36, 0.12, and 0.05 mPa, 
respectively, which should be compared with much larger mean deviations between the BS- 
BAHM approach and experiment at these separations (5, 2, 0.8, and 0.4 mPa, respectively). 
This demonstrates that the uncertainties in separations discussed in Refs. j^, Q do not 
affect the conclusion of Ref. j(| that the BSBAHM approach is excluded experimentally. 

Note that even the comparison of the static experiment in Sec. III of Ref. Q] with the 
BSBAHM approach is incorrect. Only one point at a separation z = 200 nm in Fig. 3 of 
Ref. is considered. Instead of using the original measured force vàlues, Ref. [ÏJ uses a 
maximum (not a mean) difference of -1 pN between the theoretical (as calculated in Ref. Q| 
at zero temperature) and experimental vàlues of the force. The authors of Ref. jïj explain 
this difference by the existence of a thermal correction equal to lpN. They compare this 
lpN with their predicted thermal correction of 2.56 pN at z — 200 nm and conclude that 
the result "is encouraging" . In fact, however, they have a deviation of 1.56 pN between their 
theory and the static measurement instead of a -1 pN deviation between the traditional 
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theory of Ref. M and the same measurement. Ref. [l||, however, does not inform the reader 

n n 

that in a later publication [6] the preliminary theoretical result of Ref. |10|] was recalculated 
using a more precise roughness correction (this is clearly explained in the left column on 
p.2 of Ref. Q). It was demonstrated that the static measurement of the Casimir force F^ xp 
is in fact in agreement with the theoretical vàlues given by traditional theory with a 
more precise roughness correction (see Fig. 1). Thus, if one removes this misunderstanding, 
Ref. |l| must reconcile the zero mean deviation between the traditional theory and static 
experiment (see Fig. 1) with an extra 2.56 pN thermal correction predicted by the BSBAHM 
alternative approach. This leads to the evident failure of their approach. 

Moreover, the BSBAHM theoretical approach disagrees significantly 15 
modern measurement of the Casimir force between Au surfaces of a plate and a spherical 
Iens by means of a torsion pendulum In Ref. 14 1 the experimental data were found 



131 ] with the first 



to be consistent with the theoretical Casimir force between ideal metals. A net deviation 
between the Casimir forces at a temperature T = 300 K and at a separation z = 1 /im, 
computed for ideal metals and using the BSBAHM approach, is about 25% of the Casimir 
force between ideal metals (recali that for ideal metals at z = 1 //m, T = 300 K the thermal 
correction is equal to only 1.2% of the zero-temperature force). Of this deviation, 19% is 
due to the large thermal correction predicted by BSBAHM. In spite of the fact that the 
experimental uncertainty in Ref. at 1 /im is less than 10%, the effect predicted by the 
BSBAHM approach was not observed. No mention of this important experiment is made in 
Ref. |l|. 

Quite recently the dynamic experiment of Ref. [f| was repeated (see Ref. jï^) with 
many important improvements, including a significant suppression of the surface roughness 
on the interacting surfaces, and a decrease by a factor of 1.7 (down to Az = 0.6 nm) of 
the experimental error in the measurement of the absolute separations between the zero 
roughness levels. An improvement in detection sensitivity, together with a reduction of the 
coupling between the micromachined oscillator and the environment, yielded measurements 
at smaller separations between the test bodies (160 nm instead of 260 nm). 

From the new results in Ref. ^ the conclusion was drawn that the BSBAHM approach 
to the thermal Casimir force is excluded experiment ally in the separation region from 170 nm 
to 700 nm at 95% confidence. In the separation region from 300 nm to 500 nm the BSBAHM 
approach is excluded experimentally at even higher confidence of 99% [15]. Here we illus- 



5 



trate these conclusions in Fig. 2a where the differences of the theoretical (calculated in the 
BSBAHM approach) and experimental Casimir pressures are plotted versus separation for 
14 sets of measurements containing 4066 experimental points. By contrast, in Fig. 2b the 
same differences are plotted in the case that the theoretical Casimir pressures are calculated 
using the impedance approach. In both figures the sòlid lines represent the 95% confidence 
interval for the differences between theoretical and experimental Casimir pressures as a func- 
tion of separation. It should be particularly emphasized that this confidence interval takes 
into account all experimental and theoretical errors, including in full measure the previously 
discussed error in the Casimir pressures due to experimental errors in separation distances 
as sueeested in Refs. 00. The comparison of Figs. 2a,b clearly demonstrates that the BS- 
BAHM approach is excluded by the improved dynamic experiment measuring the Casimir 
pressure, whereas the impedance approach is in excellent agreement with experiment. The 
traditional approach using the dielectric function of the plasma model is also consistent with 
the data (see Ref. [l^ for more details). 

Thus, the authors' approach is in contradiction not only with a dynamic experiment 
by means of a micromechanical oscillator and the torsion pendulum experiment ^| 
(which they leave out of their discussion), but also with a static experiment [6j and an 
improved dynamic experiment [15j which measure the Casimir force and pressure, respec- 
tively. It should be emphasized that this conclusion cannot be refuted by introducing an 
unaccounted systematic error in the measurement of the surface separation which might be 
present in the experiments of Refs. 0] in addition to the ones discussed above. The 
reason is that the influence of such an error (if it existed) decreases with an increase in 
separation whereas the contribution of the thermal correction, as predicted in the BSBAHM 
approach, increases with separation at moderate separations. Bearing in mind that the 
BSBAHM approach significantly disagrees with experiment in a wide separation range for 
several different experiments, it is easy to check that no unaccounted fixed systematic error 
is capable of bringing this approach into agreement with data within the whole range of 
measurements. As an example, in Fig. 3 we present the pressure differences (P^ 1 ' 1 — P^ xp ) 
versus separation for one typical set of measurements where all separations are decreased by 
lnm as suggested recently in Ref. |l6j]. From Fig. 3 it is clearly seen that such hypothet- 
ical systematic error is incapable of bringing the data into agreement with the BSBAHM 
approach within a wide separation range from 240 nm to 700 nm. Note also that the first 
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version of preprint lloí obtains slightly larger magnitudes of the Casimir pressure than those 
computed in Ref. |l5j] in the framework of the BSBAHM approach. (The greatest difference 



is at z = 160 nm, T = 300 K where, according to Ref. 16], P c 



ithA 



nn 



132 mPa against 



-1125. 5mPa in Ref. [l^-) The reason for this difference is that Refs. l|, [l6| used data for 
s{iÇ) from the work of Ref. 17 1. The computations in Ref. contain an error in a con- 
version coefficient from eV to rad/s (on p.313, left column the value of 1.537 x 10 15 is used 
instead of the correct value 1.51927 x 10 15 I n addition, the computations of Ref. Ï?J 

utilize Q| slightly different vàlues of the Au plasma frequency and relaxation parameter 

in Ref. (15J [u p = 9.00 eV, v(T = 300 K) = 0.0350 eV]. When both these changes are taken 



9.03 eV, u(T = 300 K) = 0.0345 eV] than were cited in Refs. and actually used 



into account in the computations of Refs. [I|,ll6j, the result P. 



c 



-1125.5 mPa of Ref . 
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is recovered. (We note that the abovementioned error does not invalidate the computations 
of the reduction factors in Ref. ^|, which were performed up to only two significant figures.) 

b) The second main claim contained in Ref. is that the approach by BSBAHM is 
consistent with thermodynamics. The calculations of are based on the Lifshitz formula 
for the free energy of a fluctuating field The dielectric permittivites of Au and Cu at 
nonzero imaginary Matsubara freqüències are taken from Ref. j^|. The contribution of the 
zero-frequency term is obtained by the substitution of the Drude dielectric function along 
the imaginary frequency axis, 

««) - 1 + wéïrTY (1) 

into the Lifshitz formula. This substitution leads to an absence of the zero-frequency con- 
tribution of the transverse elèctric mode, and, in application to metals, results in a serious 
inconsistency with thermodymanics. The point is that for perfect lattices with no defects 
or impurities but with a nonzero relaxation z/(T) and finite conductivity for nonzero T, the 
Bloch-Grüneisen law leads to v{ti) = 0. (This important property for perfect lattices remains 
even when the effects of electron-electron col·lisions are included.) As proven analytically in 
Ref. 0], for such perfect lattices the approach by BSBAHM leads to a nonzero entropy of 
a fluctuating field at zero temperature given by, 

k B 



S(z,0) 



IQ-kz 2 



ydy ln 



1 - 



y - 



JJ 



H z2 , 7( 2 

r 2 ' y 



H z2 i „2 

r 2 i y 



< 0, 



(2) 
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which depends on the parameter of the system under consideration, i.e., on the separation 
distance z (ks is the Boltzmann constant). 

In an attempt to avoid this serious problem, the authors of use a nonzero value of 
the Drude relaxation parameter at zero temperature arising from the presence of impurities. 
They can then obtain in Secs. IIA and IV a zero value of entropy at zero temperature (a 
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22|). This, however, does not solve the 



result obtained first by Bostròm and Sernelius 

inconsistency of the BSBAHM approach with thermodynamics, as is claimed in Ref. 
because it is still violated for a perfect lattice having finite conductivity and nonzero relax- 
ation at any T > 0. The unstated assumption of the authors is that perfect crystals with 
no defects or impurities do not exist and, therefore, that the laws of thermodynamics need 
not apply to them. This assumption is unphysical. Nernst and Planck formulated their 
famous theorem for the case of a perfect lattice which is truly an equilibrium system. Later 
this theorem was proven in the framework of quantum statistical physics for an y sy stem 
with a nondegenerate dynamical state of lowest energy (see, for instance, Refs. |3, 
Consequently, it is vàlid for both perfect lattices and lattices with impurities. The violation 
of the Nernst heat theorem for a perfect lattice would lead to the failure of the theory of 
electron-phonon interactions, and to the eventual abandonment of much of condensed mat- 
ter physics, statistical physics and thermodynamics. For this reason the approach advocated 
by the authors is, in fact, in violation of thermodynamics. 

In addition to the above serious problems, Ref. |lj contains several misleading statements, 
including the following: 

1) The authors discuss the so-called Modified Ideal Metal model (MIM), but do not 
mention that it violates thermodynamics. Namely, for the MIM model the free energy in 
the case of two parallel plates [Eq. (3.2) of Ref. Q]] can be identically rearranged to the form 

f MIM {z, T) = T lM {z, T) + ((3)k B T/(16Trz 2 ), (3) 

where T IM is the Casimir free energy for the usual ideal metal (2^, and ((z) is the Riemann 
zeta function. By differentiating both sides of Eq. (JHJ with respect to T one obtains 

S MlM (z, T) = S m (z, T) - C(3)/W(167r;2 2 ), (4) 

where S M1M and S lM are the entropies for the MIM and IM, respectively. Taking into account 



S 



that for the usual ideal metal S 1M — > when T — > 



we come to the conclusion 



5 M1M (a, 0) = -C(3)/V(167rz 2 ) < 0, (5) 

i.e., the MIM model violates the Nernst heat theorem. Remarkably the result for real metals 
I 

obtained in Ref. jlj coincides with that for MIM (which violates thermodynamics) at large 
separations and does not coincide with the classical limit based on Kirchhoff's law [27 1. 

2) In Sec. IV of Ref. |j| it is claimed that "a transverse elèctric zero mode ... should 
not be present according to Maxwell's equations of electromagnetism" . This is not correct. 
Maxwell's equations alone do not lead to a contribution of the transverse elèctric mode at 
zero frequency, unless they are supplemented by an adequate characterization of the mate- 
rial boundaries. More importantly, the characterization by means of the Drude dielectric 
function, used in Ref. [ïj], is inadequate to describe virtual photons. 

3) In Sec. VI the authors repeat their argument of Ref. 4J that the exact impedances, 
which depend on transverse momentum, lead to a zero contribution of the transverse elèctric 
mode to the Casimir force, as does the Drude dielectric function. Ref. Q] claims that Refs. [?], 



where the nonzero contribution to this mode was obtained, comp 



etely disregard the 



201 ] the dependence 



transverse momentum dependence. This, however, is incorrect. In Ref. 
of the impedance on a transverse momentum is considered in detail. As was demonstrated in 
Ref. [20] , this dependence disappears in the limit of zero frequency if the dispersion equation 
for the determination of photon proper freqüències is taken into account in the boundary 
conditions, and, as a consequence, the reflection properties of virtual photons on a classical 
boundary coincide with those of real photons. The transverse-momentum dependence in the 
impedance is, however,preserved at nonzero freqüències and in the reflection coemcients at 
all freqüències. Ref. |l| fails to account for this important property, thus resulting in a 
violation of the Nernst heat theorem. Ref. recalis that the Lifshitz formula contains not 
real freqüències but the i mag inary Matsubara freqüències for which there is no mass-shell 
equation. However, Ref. jl6| disregards the fact that the impedance boundary conditions 
are prior to the Lifshitz formula. They are imposed on any real or virtual electromagnètic 
wave of real frequency. The next step is the derivation of the dispersion equation and 
only after the application of the argument principle do we arrive at the Lifshitz formula 



witb the imaginary Matsubara freqüències 0. If we adniit thai the reflection properties of 
virtual photons on a real boundary coincide with those of real photons, we must impose 
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the dispersion equation on the impedance function in the boundary çondition formulated 
in terms of real freqüències. Therefore, the speculations in Ref. [l6| that for imaginary 
freqüències the angle of incidence becomes meaningless are misleading. 

4) In the end of Sec. VI of Ref. Q, the claim of Ref. |5| is repeated against the extrap- 

region to zero frequency. No 

mention is made of Ref. |28( containing the justification for this extrapolation, and demon- 



strating that the treatment of the zero-frequency mode as in Ref. [5|] results once again in a 
violation of the Nernst heat theorem. Preprint lloí claims that according to Ref. 15J real 

n 

data cannot be used in the description. Ref. 15 , however, uses exactly the same real data 

n n 

and in the same frequency range as Ref. |1J does. In fact, the disagreement between Refs. [1| 
and is only in the method of extrapolating real data to zero frequency. 

To conclude, the theoretical approach of Ref. (jj is excluded by four already performed 
experiments, namely, by the measurements of the Casimir force using a torsion pendulum 
12, 13lll4|, a micromechanical oscillator fel, and by using two determinations of the Casimir 

f n n 

pressure in a micromechanical system jy, [15J]. The main results of Ref. |l| are also in 
contradiction with fundamental physical principies such as the laws of thermodynamics. 
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FIG. 1: Differences of theoretical and experimental Casimir forces between the sphere and plate 
versus separation in the static experiment calculated and measured in Ref. (|. 
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FIG. 2: Differences of theoretical and experimental parallel plate Casimir pressures versus sepa- 



ration obtained from the improved dynamic measurement of Ref. 



P^ 1 ' 1 are calculated as in Refs. 



The theoretical vàlues for 

nn 



DDDDD (a), and for Pç as in Refs. |d,Uj 

(b). Sòlid lines 

represent the 95% confidence interval. 
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FIG. 3: Differences of theoretical and experimental parallel plate Casimir pressures versus separa- 
tion obtained from the improved dynamic measurement of Ref. after all separation distances 
have been decreased by 1 nm to take into account the hypothetical systematic error suggested by 
Ref. [Rj. The theoretical vàlues for P c ' are calculated as in Refs. Sòlid lines 

represent the 95% confidence interval. 
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